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The construction of the Agua Negra tunnels that wiU hnk Argentina and Chile under the Andes, 
the world longest mountain range, opens the possibility to build the first deep underground labo- 
ratory in the Southern Hemisphere. This laboratory has the acronym ANDES (Agua Negra Deep 
Experiment Site) and its overburden could be as large as ~ 1.7 km of rock, or 4500 mwe, providing 
an excellent low background environment to study physics of rare events like the ones induced by 
neutrinos and/or dark matter. In this paper we investigate the physics potential of a few kiloton 
size liquid scintillator detector, which could be constructed in the ANDES laboratory as one of 
its possible scientific programs. In particular, we evaluate the impact of such a detector for the 
studies of geoneutrinos and galactic supernova neutrinos assuming a fiducial volume of 3 kilotons 
as a reference size. We emphasize the complementary roles of such a detector to the ones in the 
Northern Hemisphere neutrino facilities through some advantages due to its geographical location. 

PACS numbers: 14.60.Lm,14.60.Pq,13.15.+g,95.85.Ry 



I. INTRODUCTION 

After the pioneering neutrino experiments performed 
at Homestake (South Dakota, USA) [T] and Kamioka 
(Gifu, Japan) [2^, and the great achievements by Super- 
Kamiokande [3], KamLAND [4], both in Kamioka, and 
SNO [5] (Sudbury, Canada) experiments which provided 
strong evidence of neutrino oscillation, it has been well 
recognized that deep underground laboratories can offer 
an excellent environment for neutrino experiments as well 
as for a variety of interesting scientific programs which 
include several different fields, from particle physics, as- 
trophysics, nuclear physics, to biology, geology and geo- 
physics. See [B] for a review of the world's underground 
laboratories. 

Experiments searching for very rare events, such as the 
ones induced by neutrinos or dark matter interactions, 
proton decay or performing low energy nuclear cross sec- 
tion measurements, cannot be carried out at the Earth's 
surface mainly due to the backgrounds induced by cosmic 
rays. For these experiments a reduction of the cosmo- 
genic backgrounds is crucial. This can be accomplished 
by having enough rock overburden, making clear the rea- 
son for going deep underground. 

Recently it was proposed j7j to build the first under- 
ground laboratory in the Southern Hemisphere by dig- 
ging a cave off one of the two 14 km long Agua Negra 
tunnels. These tunnels will be constructed under the An- 
des, the longest continental mountain range in the world, 
to connect Chile's region IV and Argentina's San Juan 
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province. They will provide a link between the port of 
Coquimbo, Chile (on the Pacific Ocean), and the port of 
Porto Alegre, Brazil (on the Atlantic Ocean), and several 
nearby cities, in order to facilitate trade between Asia 
and MERCOSUR (Mercado Comun del Sur or Common 
Southern Market in English) which is an economic and 
political agreement among Argentina, Brazil, Paraguay 
and Uruguay. While the exact location of the laboratory 
inside the tunnel is still under study, the rock overburden 
could be as large as ^ 1.7 km, allowing to significantly 
reduce backgrounds from cosmic ray origin. The name 
given to the proposed laboratory is ANDES (Agua Negra 
Deep Experiment Site). 

If such an underground laboratory is constructed, it 
could provide a variety of interesting scientific oppor- 
tunities for dedicated studies of neutrinos, dark mat- 
ter searches and nuclear astrophysics among other 
things. The current preliminary design of the ANDES 
laboratory is as follows [S]. There will be two large ex- 
perimental halls with dimensions of 21 x 23 x 50 m"^ and 
16 X 14 X 40 m'^, and one smaller hall for offices and mul- 
tidisciplinary experiments with size of 17 x 15 x 25 m'^. 
In addition, there will be two experimental pits, one is a 
smaller ultra-low radiation pit with 8 m of diameter and 
9 m of height, the other is a large single experiment pit 
with 30 m of diameter and 30 m of height. We are par- 
ticularly interested in this large pit where a liquid scin- 
tillator detector could be installed and used in a possible 
neutrino program for the ANDES. 

As demonstrated by KamLAND [Ij , Borexino [HI HH] , 
as well as the recent ^13 reactor experiments, Double 
Chooz [11 , Daya Bay [l| and RENO [13 , hquid scintil- 
lator detectors have a very good capability to observe 
through the inverse beta decay reaction, +P — >■ n-\-e^. 
They also can work with a low energy threshold and pro- 
vide good energy resolution. As a possible candidate for 
the ANDES neutrino detector one could consider a Kam- 
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LAND fi", Borexino [H] or SNO+ [H] like liquid scin- 
tillator detector with a fiducial mass of a few kt. In this 
paper, we assume a liquid scintillator detector based on 
alkyl benzene (C6H5C12H25), which will be used for the 
SNO+ detector, with a fiducial mass of 3 kt, containing 
^ 2.2 X 10"^^ free protons as targets, as our reference neu- 
trino detector at the ANDES, unless otherwise stated. 

We focus on the detection of neutrinos originating from 
some radioactive elements inside the Earth, the so called 
geoneutrinos, and neutrinos coming from a core collapse 
supernova (SN) in our galaxy (hereafter SN implies core 
collapse supernova). For reviews on these subjects, see, 
for example, Refs. [TH |T7] for geoneutrinos and Ref. [TH] 
for SN neutrinos. Some preliminary results of this work 
can be found in Refs. [I9H2T] . 

The first successful observations of geoneutrinos by 
KamLAND [22] and Borexino [23] open a new window 
to study the chemical composition of the Earth. It 
is believed that the geoneutrino flux has a local de- 
pendency [16) . hence having more detectors in different 
parts of the Earth capable of measuring such neutrinos 
is highly welcome. Since the ANDES laboratory is sur- 
rounded by a thick continental crust, we expect a geoneu- 
trino flux larger than at Kamioka or Gran Sasso, which 
is interesting to confirm experimentally. Moreover, com- 
pared to other existing underground laboratories, there 
are very few nuclear reactors around the ANDES location 
- a valuable advantage as reactor neutrinos are one of the 
main backgrounds for the measurement of geoneutrinos. 

After the historical observations of SN neutrinos from 
SN1987A occurred in the Large Magellanic Cloud by 
Kamiokande 2J, 1MB [211, ^^'^ Baksan [25] detectors, 
it is well understood that such SN neutrinos can play a 
very important role in uncovering the physics of super- 
novae as well as some properties of the neutrinos them- 
selves, such as mass, lifetime, magnetic moment, etc |26| . 
The low rate of nearby SN, which occurs close enough to 
the Earth so that it can be observed also by neutrinos, 
is another strong reason for having as many simultane- 
ously operating neutrino detectors as possible, in order 
to take advantage of such a rare opportunity. The addi- 
tional new neutrino detector would also help in making 
a quick alert to astronomers about the occurrence of a 
nearby SN event through the SuperNova Early Warning 
System (SNEWS) network [27|. 

The ANDES neutrino detector, if constructed, can cer- 
tainly make some relevant contribution to SN neutrino 
observations. Furthermore, the location of the ANDES 
laboratory in the Southern Hemisphere can provide a bet- 
ter chance to observe the Earth matter effect for SN neu- 
trinos by combining with other detectors in the Northern 
Hemisphere. If Earth matter effect is observed, the neu- 
trino mass hierarchy may be determined and, at the same 
time, evidence that different SN neutrino flavors manifest 
significantly different energy spectra can be provided. 

The organization of this paper is as follows. In Sec. 
jnlwe discuss in detail the potential of the ANDES neu- 



trino detector for geoneutrino observations. In Sec. Ill we 



discuss SN neutrino observation at the ANDES neutrino 
detector. Sec. IIVI is devoted to discussions and conclu- 
sions of our results. While we follow previous works for 
most of the calculations done in this work, for the sake of 
completeness and to be self-contained, we describe some 
details of our numerical calculations in the Appendices. 



II. GEONEUTRINOS 

A. Introduction 

The deep interior of the Earth, governed by high pres- 
sure and temperature, is the last frontier on our planet, 
which has not yet been explored by a human being. The 
deepest hole ever made so far is 12.3 km down from the 
Earth's surface [16], only about 0.1% of the Earth's diam- 
eter, and so even the top of the mantle has not yet been 
reached. In the near future, a direct access to the upper 
mantle will be possible as part of the missions of the in- 
ternational scientific research program named lODP (In- 
tegrated Ocean Drilling Program) [55] . However, to reach 
the lower part of the Earth's mantle, located at a depth 
of about 660 km from the surface, seems to be an impos- 
sible task. 

So far, there are basically two different approaches to 
overcome the direct inaccessibility of the Earth's under- 
ground below 10 km: seismology and geochemistry. Seis- 
mological data permit us to indirectly reconstruct the 
matter density profile of the whole Earth. Geochemistry, 
however, can only access the composition of the Earth 
close to the surface. For that one uses various rock sam- 
ples coming from the Earth's crust as well as very limited 
samples from the top of the mantle, thanks to volcanic 
activity and orogeny. 

In 2005 the KamLAND experiment [l^ reported for 
the first time the detection of Dg coming from the de- 
cay chains of the radioactive isotopes 238y ^^^^ '^^'^Th 
in the Earth, by using the inverse beta decay reaction 
+ j3 — 7- n + e+ . These geoneutrinos can provide use- 
ful information not only relevant to the Earth's interior 
chemical composition but also shed light on the source 
of the terrestrial heat production, opening a window to a 
new scientific field, neutrino geophysics or neutrino geo- 
science. In 2010, another experiment, Borexino, located 
in the Gran Sasso Laboratory in Italy also reported the 
measurement of geoneutrinos [23j , further contributing to 
the start of neutrino geoscience. A partial list of previous 
works on geoneutrinos can be found in Refs. [THl [2M55] . 

It is estimated that the entire Earth generates about 
40 TW, corresponding to ~ 10,000 reactors. It has been 
considered that most (or all) of the heat is generated by 
the energy deposited by the decay of radioactive elements 
like U, Th and K in the Earth's interior. By measuring 
geoneutrinos, which are the direct product of such de- 
cays, one can infer the total amount of U and Th inside 
the Earth. We note that the amount of K can not be 
inferred directly by the current detection method since 
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FIG. 1: Isocontour map of the Earth's crust thickness (in km) in the plane of 9 (latitude)-(^ (longitude) based on the model 
found in Ref. |35] . The positions of the Gran Sasso, Kamioka, Sudbury and ANDES underground laboratories are also indicated 
by asterisks. 



the energy of geoneutrinos coming from K is below the 
threshold of the inverse beta decay reaction. 

The measured geoneutrino flux, reported recently by 
the KamLAND experiment [31], is 4.3;i[;^ x 10^ cm-^g-i. 
By taking into account neutrino oscillations, this corre- 
sponds to a total emitted flux of g x 10^ cm~^s~-^. 
On the other hand, Borexino experiment [10] reported 
in terms of the observed number of events, 3.9^]^'3(l3'2) 
evcnts/(100 ton yr) at 68% (99.73 %) CL. This result is 
about 70 % higher than that obtained by KamLAND [33] , 
though both results are consistent with each other within 
the experimental and theoretical uncertainties. 

By combining the results from KamLAND and Borex- 
ino, the observed geoneutrino flux corresponds to a heat 
production of 20l:g g TW [33]. While Borexino seems to 
favor, KamLAND results disfavor the so called fully ra- 
diogenic model, the model where all the terrestrial heat 
comes from the decay of the radioactive elements in the 
Earth crust and mantle. KamLAND alone disfavors this 
model at 98.1% CL while the combined data of these 
two experiments slightly reduces the significance of this 
rejection to 97.2% CL [34] . 



B. Why measure geoneutrinos in the ANDES? 

For various reasons it would be very interesting if the 
measurement of geoneutrinos can also be done at the 
ANDES Laboratory. First we note that the location of 
the ANDES laboratory, 30° 15' S and 69°53' W, is sur- 
rounded by the Andes mountain range which means that 



the thickness of the crust around the laboratory is signif- 
icantly larger than the average Earth crustal thickness, 
leading to an expected larger geoneutrino flux. This is be- 
cause in the Earth's crust the concentration of U and Th 
is expected to be significantly larger than in the deeper 
mantle. 

In Fig. [T] we present the isocontour map of the Earth 
crust thickness based on the model found in Ref. (35) . 
In Fig. [2] we present the magnified version of Fig. [l] 
around the ANDES laboratory. From these figures, we 
can appreciate the difference of the local crust thickness 
around the ANDES laboratory in comparison with other 
locations on our planet. Indeed, roughly speaking, the 
expected geoneutrino flux at the ANDES laboratory site 
is larger than that for KamLAND and Borexino by about 
30 % and 20 %, respectively, as we will see below. It is im- 
portant to confirm such local dependence of the geoneu- 
trino fluxes. Because of such a local (site) dependence, 
it would be also very interesting to measure the geoneu- 
trino flux at a location surrounded by oceanic crust, such 
as in Hawaii |36j . where the main contribution to the 
geoneutrino flux is expected to come from the decay of 
radioactive elements in the mantle. 

Second, there are very few nuclear reactors near the 
ANDES laboratory which is a great advantage for mea- 
suring geoneutrinos. As it is well know from the re- 
sults of the KamLAND experiment, low energy pro- 
duced by nuclear reactors is one of the most impor- 
tant backgrounds for geoneutrino observation. Indeed, 
the Borexino detector, despite being smaller than Kam- 
LAND, demonstrated so far a better performance than 
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FIG. 2: Same as in Fig.[l]but around the ANDES laboratory. 

KamLAND as far as geoneutrino measurements are con- 
cerned. This seems to be due to the fact that Borexino, 
located in the middle of Italy where there are no nearby 
reactors, is exposed to a much lower Df, background from 
reactor origin. 

Taking into account the nearest Argentinian nuclear 



reactors, the Embase 2.1 GW (thermal power) as well as 
the Atucha I 1.2 GW and Atucha II 2.1 GW reactors, 
which are located, respectively, 560 km and 1080 km 
away from the ANDES Laboratory, we have estimated 
a total reactor background of 8.8 events/yr/3 kt and 2.2 
events/yr/3 kt in the geoneutrino energy range. These 
numbers are given in the absence of neutrino oscillations, 
which will further reduce them if taken into account. 
Though we have considered only the contribution from 
these nearby reactors, our estimation is similar to the one 
which can be inferred from Fig. 2 of Ref. [37] which shows 
the world map of the isocontours of the expected number 
of events induced by neutrinos coming from 201 nuclear 
reactor power stations all over the world. The reactor 
neutrino background we found for ANDES is more than 
10 times smaller than the one expected for Borexino, 5.7 
events/yr/100 tons (total number of reactor induced 
events in the presence of oscillation) [TOj . 

C. Calculating the geoneutrino flux 

We follow our (two of the authors of this paper) pre- 
vious work |33j with some update and improvements, in 
order to compute the differential flux of 9^ produced in 
the decay chain of radioactive isotopes ^sSjj ^^^^ '^^'^Th 
that will be measured at a detector position r on the 
Earth, which can be expressed by the following integral 
performed over the Earth's volume V^, 



E 

fe=U,Th 
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p(r') Ck{r')nk 
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r-'|2 



Tk-rrik 



X Pp^(^p^,|r-r'|) X fk{E^J, 



(1) 



where p(r) is the matter density, Cfc(r), r^, and Uk 
are, respectively, the mass abundance, half-life, atomic 
mass and the number of emitted per decay chain cor- 
responding to element k =^^^U, ^^^Th, fk{EpJ is the 
normalized spectral function for element k [38] . 

We will assume the concentrations of U and Th, Cfc(r), 
take different values in the Earth's crust and mantle lay- 
ers as given m Table [!} These are our reference values 
which are based on Ref. [17]. Oceanic and continental 
Earth's crust are divided, respectively, into two and four 
layers whereas the Earth's mantle is divided into two lay- 



ers (see Table |T]). We assume no U and Th in the Earth's 
core. The so called fully radiogenic model assumes higher 
U, Th and K abundances in the mantle in order that the 
total observed Earth's heat can be fully explained by the 
decay of these radio active elements. We have used the 
Earth crust model taken from Ref. [35] (shown in Figs, [l] 
and[2|. 

{Ei>^ , |r— r'l) describes the survival probability of De 
produced at r' but measured at r, which can be averaged 
out, as a good approximation, and bring out from the 
integral the term: 



{PcJ ~ (sin'* 9i3 + cos"* ^13 1 - sin^ 29i2 sin 
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iE. 



r r 



sm f^is-fcos 



713 



1 ~ - sin^ 2^12 ~ 0.55, (2) 



r 



where Am^i = — rn^ 



being the neutrino masses, sin 9i2 



7.5 X 10"^ eV"^, mi, TO2 
2^ - - 0.31 and sm^Oi 



'13 



0.025 [39]. In this work we use the standard neutrino 
mixing parameterization found in Ref. [40] . 
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Crust Crust 
Mantle Mantle 



Layer 


cu (/i g/g) 


CTh (m g/g) 


Oceanic Sediment 


1.68 


6.91 


Oceanic Crust 


0.1 


0.22 


Continental Sediment 


2.8 


10.7 


Upper Continental Crust 


2.8 


10.7 


Middle Continental Crust 


1.6 


6.1 


Lower Continental Crust 


0.2 


1.2 


Upper Mantle 


0.012 


0.048 


Lower Mantle 


0.012 


0.048 



TABLE I: U and Th mass abundances in different layers of 
the Earth's crust and mantle used in this work [IT]. In the 
Earth's core these mass abundances are assumed to be zero. 

Cumulative Flux at the ANDES 
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FIG. 3: Cumulative flux of geoneutrinos coming from the 
decay of ^^^U as a function of the distance from the ANDES 
neutrino detector. 



We have computed the expected total fluxes of at 
ANDES coming, respectively, from U and Th to be. 



$Th 



-1 (3.04 X 10^ 



cm 



-1) and 



4.78 X 10^ cni~2 g-i (2.6O x 10^ cm^^ g-i) with- 
out neutrino oscillation (with oscillation) . In Fig. |3] we 
present the geoneutrino cumulative flux for the U chain 
as a function of the distance from the ANDES labora- 
tory. We observe that 50% of the flux comes from ~ 200 
km from the detector and about 20% of the flux comes 
from the mantle. 

In Fig. |4] we show our expectations for the total os- 
cillated geoneutrino flux at Kamioka, Gran Sasso, SNO, 
Hawaii, Pyhasalmi and ANDES, discriminating the crust 
and mantle contributions in each case. Pyhasalmi in Fin- 
land is a possible site for the proposed 50 kt LENA neu- 
trino detector [51]. We also show KamLAND [33] and 
Borexino data [231 points to compare with the precision 




Kamioka Gran Sasso SNO Hawaii Pyhasalmi ANDES 

FIG. 4; Total geoneutrino flux (oscillated) expected at 
Kamioka, Gran Sasso, SNO, Hawaii, Pyhasalmi and ANDES. 
We show the expected contribution from crust and mantle, in 
each well as KamLAND 34 and Borexino data [23] 

points. We also show (in blue) the precision of the expected 
measurement by ANDES after 5 years at 1 and 3 a CL. 



of the expected measurement by ANDES after 5 years 
of data taking. According to Ref. [32] the current Kam- 
LAND and Borexino results combined imply the geoneu- 
trinos from the mantle have been observed at 2.4 a CL. 
Clearly ANDES by itself, after 5 years, is able to estab- 
lish the mantle geoneutrino component at a level of about 
3 cr or better. 

Let us now discuss the expected number of geoneutrino 
induced events at the ANDES neutrino detector. For one 
year of operation (3 x 10'' s) and 80% detector efflciency, 
we have calculated the total number of geonetrinos ex- 
pected at the ANDES reference detector to be 82.4 (64.8 
from U, 17.6 from Th). About 16 of these events would 
be from the mantle and 35 events would have En > 2.3 
MeV, coming exclusively from the U chain. To illustrate 
the site dependence we show in Table jll] our estimation 
for the corresponding number of geoneutrino events in 
different locations assuming the same reference detector 
(the same number of free protons, efficiency and expo- 
sure). From this table we see that the expected number 
of geoneutrino events at the ANDES location is compa- 
rable to SNO and to Pyhasalmi. 

Such a detector operating during 10 years could accu- 
mulate more than 800 geonetrino events (160 from the 
mantle alone), allowing not only for a better determi- 
nation of U and Th mass abundances in the crust and 
mantle but also for the investigation of their presence in 
the Earth's core. Clearly if an even larger detector, say 
10 kt, could be envisaged the scientific reach could be 
even more significant. 
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Location 


Number from U 


Number from Th 


Total 


Gran Sasso 


53.8 


14.7 


68.5 


Kamioka 


45.7 


12.4 


58.1 


Hawaii 


18.5 


5.0 


23.5 


Sudbury 


63.2 


17.2 


80.4 


Pyhasalmi 


66.1 


18.0 


84.1 


ANDES 


64.8 


17.6 


82.4 



TABLE II; Expected number of geoneutrino events for our 
reference 3 kt liquid scintillator detector operating during a 
year with 80% efficiency at different locations. 



III. SUPERNOVA NEUTRINOS 
A. Preliminaries 

If we consider the entire universe, the SN event rate is 
not so low, several SN explosions per second. However, if 
we restrict to our galaxy, more interesting in terms of SN 
neutrino observations, the estimated event rate of such 
nearby SN drops down to about < 3 per century 
|45] . In fact, in last ^ 30 years, since when the Baksan 
neutrino detector started to operate in 1980, only neu- 
trinos from the explosion of SN1987A in the Large Mag- 
ellanic Cloud, one of the satellite galaxies of the Milky 
Way, were observed by the Kamiokande [2], 1MB [24], 
and Baksan [5S] neutrino detectors. 

Fortunately today, compared to that epoch, we are 
better prepared for SN neutrino observations as there 
are much larger neutrino detectors, such as Super- 
Kamiokande [15] and IceCube 07], which are currently 
in operation. However, since the nearby SN rate is quite 
low, it is better to have as many neutrino detectors as 
possible, to be ready for the next SN event. This will 
maximize the chance to obtain as much information as 
possible on SN neutrinos, leading to a better understand- 
ing of the SN explosion dynamics. 

Since the last galactic SN, SN1604, was observed more 
than 400 years ago and we do not know when exactly 
the next one will occur, it is important to have neutrino 



detectors with larger running time. Having many neu- 
trino detectors is also helpful in forming a network like 
SNEWS [57] that will readily alert the astronomers about 
the occurrence of a nearby SN event enabling them not 
to miss the initial phase of the time evolution of the SN 
light curve. 

It is theoretically expected [IS] that almost all 99 
%) the energy released by gravitational collapse is car- 
ried away in the form of neutrinos, which was consistent 
with the observed data of SN1987A neutrinos. Roughly 
speaking, the neutrino emission from a SN explosion can 
be divided into four periods: (i) the infall phase, which 
starts several tens of ms before the bounce; (ii) the shock 
breakout neutronization burst, which lasts up to a few 
tens of ms after the bounce; (iii) the accretion phase, 
during from a few tens of ms up to several hundreds of 
ms after the bounce and (iv) the Kelving-Hclmholtz cool- 
ing phase, during up to ^ 10-20 s after the bounce. 

Emission of starts during the infall phase, though 
the luminosity is not yet so large. In the neutronization 
burst phase, there is a strong I'e burst in a very short 
period of time, ^ 10 ms, and the emission of the other 
flavors, as well as of Vg, is suppressed. During the ac- 
cretion phase, the fluxes of Vg and Pe are expected to 
be significantly larger than that of the other flavors and 
the energy hierarchy {E^J < (Egy^) < (E^^) is expected. 
Here we use the notation ly^ to refer to any non-electron 
neutrino since for our purpose they can be treated, in 
good approximation, as a single species. During the cool- 
ing phase the emission of all flavors of neutrinos and an- 
tineutrinos, with similar luminosities, is predicted. 

It is considered that the energy spectra of SN neutri- 
nos can be approximately described by Fermi-Dirac dis- 
tributions with some non-zero chemical potential, which 
is in general necessary to account for the non-thermal 
feature of the SN neutrino spectra. In this work, for 
any flavor, Va, we will use the following parameteriza- 
tion, which is based on the numerical simulations by the 
Garching group [iSHST] . for the SN neutrino spectra at 
the Earth in the absence of neutrino oscillation. 



E 

I 



exp 



E 



(3) 



where D is the distance to the SN, is the total num- 
ber of Va emitted, (E^^) is the average energy of and 
Pa is a parameter which describes the deviation from a 
thermal spectrum (pinching effect) that can be taken to 
be 2 — 4, r(/3a) is the gamma function. 

This parameterization seems to describe better the SN 
neutrino spectra obtained by numerical simulations. We 



note, however, that our results would not change much 
even if we had used instead Fermi-Dirac distributions 
with a non-zero chemical potential. During the neutrino 
emission, as many SN simulations indicate, the shape of 
F^^ (E) is expected to change in time, which means that 
the average neutrino energies as well as their luminosities 
are, in general, functions of time. We do not explore this 
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feature in this work ^. 

For the sake of discussion, unless otherwise expHcitly 
stated, we use the SN parameters summarized in Table 
mil as our reference values. We assume that the total 
energy released by SN neutrinos is 3x10^'^ erg, which 
is equally divided by 6 species of active neutrinos, and 
consider 10 kpc as a typical distance to the SN. The 
chance that the actual distance could be even smaller, say 
5 kpc, seems, however, to be not so small. This can be 
seen in Fig.[l6]in the AppendixjCj where the expected SN 
distributions as a function of the distance from the Earth, 
based on the distribution models considered in |52) . are 
shown. 



Reference SN parameters 
Distance to the SN : D = 10 kpc 
Spectra parmetrization in Eq.([3]l with /3a = 4 for all flavors 
{E,J = 12 MeV {Ei,J = 15 MeV {E^J = 18 MeV 
ETJ^ = {E^J$^^ = 5 X lO'^^ erg for all flavor 



TABLE III: Reference SN parameters to be used throughout 
this paper, unless otherwise explicitly stated. 

Due to neutrino oscillations in the SN envelope, the 
SN neutrino flux spectra at Earth are in general expected 
not to be given by Eq. ([3| but by flavor mixtures (super- 
positions) of them. Without loss of generality, the flux 
spectrum of observable i>e at the Earth can be expressed 
as 

F,^ (E) = piE) Fl (E) + [1 - pmFl (E) , (4) 

where F^ (E) and (E) are, respectively, the original 
spectra of Vf, and P^.^ at the SN neutrinosphere in the 
absence of any oscillation effect. Eq. Q implies that as 
long as the final observable spectrum at Earth, Fp^ (E) , is 
concerned, thanks to the very similar spectra of D^^ and D^- 
in the SN, we can work within the effective 2 flavor mixing 
scheme. Note that if there is some difference between i/^ 
and Pr at the origin, (E) must be given by c^F^^ + 
CrFp (E) where the values of and Ct, which satisfy + 
Ct- = 1, depend on the mixing and oscillation scenario. 

The function p{E) is the survival probability. This 
is, in general, a function of the neutrino energy E and 
incorporates all oscillation effects inside the SN, includ- 
ing the standard Mikheyev-Smirnov-Wolfenstein (MSW) 
effect [S3], the so called collective effects [51] as well as 
other possible effects coming from shock wave [SS], den- 
sity fluctuations [SS], turbulence [57], etc. Note thatp(F) 
can also depend on time. 



Since the observation of SN events in the ANDES laboratory 
is hkely to be statistically limited, as we will see in the next 
subsections, and also due to the SN model uncertainty on the 
time dependence of F^^{E), a detailed time dependent analysis 
will not be considered here. 



From this expression, one can see that for a given value 
of p{E) ^ 1, a larger difference between Fp^{E) and 
F^^{E) corresponds to a larger observable effect of SN 
neutrinos at the Earth. Clearly, there is no observable 
effect ifFO(F)=FOjF). 

According to recent SN simulations [SBj-jBD], especially 
during the cooling phase, the mean energies of the differ- 
ent flavors and their luminosities are likely to be similar. 
If so, as mentioned above, any observable oscillation ef- 
fect tends to vanish. Nevertheless, we should keep in 
mind that the SN neutrino spectra must be ultimately 
determined by the future SN neutrino observations, inde- 
pendently from any theoretical predictions of numerical 
SN simulations. 

If only the standard MSW effect in the SN envelope is 
operative, then it is straightforward to calculate the value 
of p{E) for a given value of ^13 and a fixed mass hierar- 
chy |53j . Since we know now the value of ^13 rather well 
thanks to the recent measurements by accelerator [651166] 
and reactor experiments [TTHT3] (see also the combined 
analysis in Ref. |67|), the only open question in this sce- 
nario is the neutrino mass hierarchy. This will be referred 
to as the standard scenario. 

The presence of other effects beyond the standard sce- 
nario, such as collective oscillations, shock waves, turbu- 
lence, etc., can cause a significant modification to p{E), in 
the case of the inverted mass hierarchy. However, some 
recent studies [B2J |S3] indicate that as long as the ac- 
cretion phase is concerned, collective effects should be 
strongly suppressed due to the high matter density. This 
implies that in this phase only the standard scenario, the 
usual MSW effect, would be operative. 

For the normal mass hierarchy, in the standard sce- 
nario, neglecting corrections from other effects, we expect 
that p{E) « [ST, and therefore, 

F,^iE)^cl,Fl{E) + sl,FliE). (5) 

On the other hand, for the inverted mass hierarchy, by 
taking into account the recent observation of non-zero 
6*13 by accelerator [65] IM] and reactor [TTJ [12] experi- 
ments, it is expected that Fp^{E) « F^^{E) [S3], i.e., 
p{E) — >■ 0, due to the adiabatic conversion inside the SN 
driven by the mass squared difference Am|2, relevant to 
atmospheric neutrinos. Hence, it is expected that inde- 
pendently of the mass hierarchy, even after taking into 
account possible corrections due to collective oscillations, 
shock waves, density fluctuations, etc., the value of p{E) 
in Eq. Q satisfy < p{E) < c?2 - 0.69 "64]. 

In Fig. [5] we show the expected neutrino flux spectra 
at the Earth for a SN located at 10 kpc from the Earth, 
using the typical average SN neutrino energies we con- 
sider in this work, (F^J = 15 MeV and (F^^J = 18 MeV 
and 22 MeV for the upper and lower panels, respectively. 
The expected De flux spectra at the Earth are shown by 
the solid red curves for the case of normal mass hierarchy 
which is given by Eq. ([5]), whereas for the inverted mass 
hierarchy, spectra is given by the dashed green curves. 
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FIG. 5: The expected neutrino flux spectra at the Earth for 
a SN located at 10 kpc from the Earth. To calculate SN 
neutrino flux at the Earth, we assumed that at the SN neu- 
trinosphere, {Ea^} = 15 MeV (dotted blue curves) whereas 
— 18 MeV and 22 MeV for the upper and lower panel, 
respectively (dashed green curves). For both cases, we show 
the expected Hux at the Earth by the solid red curves, 
given by Eq. ([5|. For simplicity, we consider that only the 
standard MSW effect plays a role in oscillations so that p = 
0.69 (0) corresponds to the normal (inverted) mass hierarchy. 



We note that since we assumed that the total SN neu- 
trino luminosity was equally divided into the 6 species of 
neutrinos, a larger average energy implies a smaller flux, 
as we can see from the dashed green curves in Fig. |5] 



B. Inverse beta decay reaction 

In this paper, we focus on two main reactions. One 
is the inverse beta decay i/g + p — > n + e+, which is the 
predominant channel due to its larger cross section, and 
the other is the proton-neutrino elastic scattering, v + 
p — + p, which is useful to determine the original 
spectra [68|, l69| . Based on the fluxes shown in Fig. [sj we 
present in Table |IV| the expected number of events we 
have computed for the ANDES neutrino detector in the 
absence and presence of neutrino oscillations, for normal 
and inverted mass hierarchy and different {Ei,J). 

Only for Table |IV[ for the purpose of comparison, we 
have considered three different chemical compositions for 
the liquid scintilator. These are mixtures that are al- 
ready used by the existing or planned detectors, Kam- 
LAND |4j, Borexino [H] and SNO-h [IS]. KamLAND 
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FIG. 6: Expected energy distribution of events coming 
from a SN located at 10 kpc from the Earth at the ANDES 3 
kt detector. We assumed that at the SN neutrinosphere, {E^^) 
— 15 MeV (dotted blue curves) whereas {E^_^) = 18 MeV and 
22 MeV for the upper and lower panels, respectively (dashed 
green curves). For both cases, we show the expected Ve Hux 
at the Earth for the normal (inverted) mass hierarchy by the 
solid (dashed) red curves. 



is based on the mixture of 80% of C12H26 and 20% of 
C9H12, whereas Borexino and SNO-j- are based on CgHi2 
(pseudocumene) and C6H5C12H25 (alkyl benzene), re- 
spectively. For the rest of the paper we assume the same 
composition of SNO-f . 

Regarding the oscillation probabilities, it is sufficient 
at this point to consider the two extreme cases of the 
standard scenario where p = 0.69 and p = Q correspond- 
ing, respectively, to the normal and inverted mass hier- 
archy. As mentioned before, apart from some possible 
non-trivial energy dependence, one expects the actual 
number of events to lay between these two cases even 
if other oscillation effects come into play. 

In Fig.|6]we show the expected v,, event number distri- 
bution dN/dE for our reference ANDES neutrino detec- 
tor in the absence and presence of neutrino oscillations. 
Depending on the oscillation probabilities, we expect to 
have ~ 800 - 1000 events for 3 kt (see Table |lV]). In the 
presence of oscillation, the energy spectrum of gets 
harder whereas its total flux will decrease, as we assume 
that the original value of (i?^^ ) is constant for all 6 
species. Nevertheless, the oscillation effect makes the ex- 
pected observed event number larger, as the cross section 
depends on ^ E'^ , which overcomes the reduction of the 
ffux due to oscillations. 
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Chemical Composition of the Scintillator 




Reaction 


(a) C12H26 + C9H12 
( 80% + 20% ) 


(b) C9H12 
pseudocumene 


(C) C6H5C12H25 

alkyl benzene 


Assumptions 


Ue+P ^ n + e'^ 


873 


630 


762 


No Oscillation 


Ue+p ^ n + e'^ 
Ve+p^ n + 


924 
1038 


669 
750 


804 
903 


p = c?2 = 0.69 (NH), (£'.J= 18 MeV 
p = 0.0(IH), (£.J=18MeV 


i>e+p ^ n + e'^ 
+ p ^ n + 


957 
1140 


690 
825 


834 
993 


p = c?2 = 0.69 (NH), {E^^}= 20 MeV 
p = 0.0(IH), (£„J=20MeV 


Ve + p ^ n + e'^ 
Ue+p ^ n + e'^ 


987 
1239 


714 
894 


858 
1080 


p = c?2 = 0.69 (NH), {E,J= 22 MeV 
p = 0.0(IH), {£,J=22MeV 


u + p ^ u + p 


294 


318 


453 


all flavors Tquo > 0.2 MeV, {E^J= 18 MeV 


u + p ^ u + p 


399 


405 


561 


all flavors Tque > 0.2 MeV, (£'^,)= 20 MeV 


u + p ^ u + p 


510 


492 


663 


all flavors Tquc > 0.2 MeV, (-B„,)= 22 MeV 



TABLE IV: Expected number of SN neutrino induced events for the inverse beta decay and proton-neutrino elastic scattering 
for the 3 types of liquid scintillators with the fiducial mass of 3 kt for a SN at 10 kpc from the Earth. They are (a) 80% of 
C12H26 and 20% of C9H12 used for KamLAND, (b) C9H12 (pseudocumene) used for Borexino, and (c) C6H5C12H25 (alkyl 
benzene) to be used for SNO-I-. NH and IH indicate the normal and inverted mass hierarchies, respectively. For iz + p — > v + p, 
we considered the kinetic (quenched) energy of the recoil proton, Tque, larger than 0.2 MeV following Refs. [68ll69] . 



Due to neutrino oscillations the observed spectra at 
the Earth is, in general, a mixture of the original ones it 
would be desirable to be able to reconstruct the original 
spectra from data. 

In principle, by fitting the observed positron spectrum 
of i/e induced events at the detector, one can try to recon- 
struct the original spectra (luminosites and average en- 
ergies) of De and Dx as done, for example, in Refs. [701171] 
but in practice, this may not be so trivial to do, even for 
a much larger detector such as Hyper-Kamiokande [75], 
due to the possible presence of degeneracies of the SN pa- 
rameters [7S] or some unexpected large deviation of SN 
neutrino spectra from what is usually assumed. 

As far as the reconstruction of the spectra of 9^ is con- 
cerned, the use of the proton-neutrino elastic scattering 
seems to be more promising [68] [69] . This will be dis- 
cussed in the next subsection. We propose to combine 
both, the inverse beta decay and proton-neutrino reac- 
tions, in order to identify the oscillation effect, or deter- 
mine SN parameters if oscillation effect is known, as we 
will discuss in subsection IIIIDl 

C. Proton- neutrino elastic scattering 

In this section, we focus on the proton-neutrino elas- 
tic scattering discussed in Refs. [H] |HH1 HH]- Since the 
proton-neutrino elastic scattering, v + p ^ + P, for all 
flavors, occurs via neutral current interactions, one can 
measure the total SN neutrino flux above a certain energy 
threshold, without worrying about any oscillation effect 
among active flavors, in a similar manner as the SNO ex- 
periment was able to measure the total solar neutrino flux 
of active flavors [5] . If the average energies of fe and De 
are significantly lower than that of v^, then by counting 
the total number of events above a certain recoil proton 



energy, one can measure the total flux (luminosity) of Vx 
flavor with reasonably good precision [SS] [5^ • If the en- 
ergy spectra of SN i^g, D^, and Vx are rather similar, one 
can then try to determine the total neutrino flux. 

To compute the number of events induced by the ly-p 
elastic scattering, we follow the analysis procedure de- 
scribed in Refs. [55] HH]- We have checked that by using 
the same information and assumption used in these ref- 
erences we could obtain results which are in good agree- 
ment with the ones presented in these works. 

In Fig. [7] we show the expected event number distri- 
bution, diV/dTquc, as a function of the quenched kinetic 
proton energy, Tquoi for the ANDES reference neutrino 
detector and for the reference SN parameters summarized 
in Table Iml 

Following [5^, the energy range of Tquo > 0.2 MeV 
was considered mainly because of the backgrounds com- 
ing from radioactive decays in the scintillator and sur- 
roundings, such as the one coming from the (3 decay of 
i-^C. The qualitative behavior of dN/dT^^^ for the AN- 
DES detector is very similar to the ones obtained in [21] 
for KamLAND, Borexino and SNO-I- detectors (see Figs. 
1-3 of this reference), the main difference is the total 
number of events. 

We have studied how accurately one can reconstruct 
the original neutrino flux by the energy distribution of 
recoil protons in the ANDES detector. As demonstrated 
explicitly in Ref. [69. by considering the quenched proton 
recoil energy larger than 0.2 MeV, one can reconstruct 
the original neutrino flux above > 25 MeV. The idea 
is to make an inversion of the calculation in order to 
produce the curves shown in Fig. [7] 

In Fig.[8]we show our result. For our reference SN, with 
the 3 kt ANDES detector, one can try to determine the 
original neutrino flux with the precision of ~ 15 % for the 
neutrino energy ^ 20-40 MeV. Note that this result does 
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Quenched Kinetic Energy T [MeV] 



FIG. 7: Distribution of events as a function of the quenched 
kinetic energy of the proton for the ANDES reference detector 
using the reference SN parameters defined in Table [Till 

not depend on the uncertainty on the effects of neutrino 
oscillations among active flavors, which could occur in the 
SN envelope. While the expected precision is not as good 
as the one expected for the proposed much larger LENA 
detector [41], which has 50 kt fiducial volume, as we can 
compare our results with the one shown in Fig. 9 of 
Ref. [BS] , the expected precision for the ANDES detector 
is better than the currently existing (planned) detectors 
like KamLAND and Borexino (SNO+). 




ho 30 40 50 60 70 80 



Neutrino Energy E [MeV] 

FIG. 8: Reconstruction of the original SN neutrino flux (in- 
dicated by the red solid circles with error bar of 1 a), from 
the simulated future data of proton-neutrino elastic scatter- 
ing in the ANDES detector under the same assumptions used 
in Fig. [T] for the reference SN parameters defined in Table 

nm 

Next let us discuss with which precision we can try to 
determine the original mean energy. For this purpose, 
we perform a analysis by considering the input (true) 



values of SN parameters for {E^J = 15, 18 and 22 MeV 
and Elf = 5 x 10^^ ^^.^ -y^g (^jg^^e E^f such that 4£'*°* 
gives the total energy carried away by non-electron flavor 
neutrinos. For simplicity, we fix the other SN parame- 
ters for and to our reference values, and only vary 
and (Ei,^) in our fit, in order to have some feeling 
about the precision we can achieve. We show our re- 
sult in Fig. [9] where the allowed parameter regions in the 
plane of and (E'^^ ) is shown. Here only the statisti- 
cal uncertainties were taken into account. From this, we 
conclude that we can determine (i?^^ ) and with the 
precisions of - 6% (25%) and - 20%" (40%), respectively, 
at 1 (3) (7 CL, for the case where the true values of (E^^) 
and were assumed to be 18 MeV and 5 x 10^^ erg, 
respectively. 
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^2 14 16 18 20 22 24 
<E^ > [MeV] 

FIG. 9: Sensitivity of the determination of the original iS^"' 
and (Ei,^) at 1,2 and 3 cr CL, by proton-neutrino scattering 
events measured at the ANDES detector. Only the statisti- 
cal uncertainties were taken into account. We consider three 
different input values for the average energy, (E^^) — 15 
MeV, 18 MeV and 21 MeV while that for Elf was fixed to be 
5 X 10^^ erg (they are indicated by the asterisk). The other 
SN parameters are defined in Table [Till 



D. Comparison of CC and NC induced events 

As mentioned in the previous subsection, one of the 
important features of the reaction induced by proton neu- 
trino elastic scattering is that it does not depend on the 
neutrino oscillation among active flavors, as this reaction 
is induced by neutral current (NC). On the other hand, 
the inverse beta decay is induced by charged current (CC) 
interactions and so does depend on the neutrino oscilla- 
tion probability, as can be seen in Table |IV[ 

Therefore, by comparing these two kinds of events, one 
can try to infer, to some extent, the oscillation effect 
or the value of p in Eq.Q in a less model dependent 
way. Or conversely, if the mass hierarchy is known (from 
some other source, like one of the long-baseline neutrino 
oscillation experiments) by the time the next galactic SN 
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occurs, and p can be predicted in advance, one can try 
to determine better some other characteristic of the SN 
neutrinos. We observe that the comparison of CC and 
NC events was also considered in [73]. 

For the purpose of illustration, let us define the follow- 
ing ratio. 



RiN^^p/N^p) 



^ UT) 



obs 



(6) 



where {Nj^^p/N^p)°^^ means the ratio of the observed to- 
tal number of events induced by the inverse beta decay 
and proton-neutrino elastic scattering. 

In Fig. [To] we show this quantity as a function of the 
true value of {E^^ ) . For p, we considered the two extreme 
values of the standard scenario, p = and 0.69. For this 
study, we consider the two cases where D = 5 and 10 kpc, 
indicated by the darker and lighter color, respectively. 

We note that if (-E^^) is significantly different from 
(-Ep^), one can try to identify the presence of oscillation 
effect by inferring the value of p with some precision, 
which could be possible especially for the case where the 
true value of (E^^) turns out to be quite different from 
{Ejy^). This may be the case during the accretion phase 
as energies and luminosities of Pe and may be signif- 
icantly different, as recent SN simulations indicate [58l - 
150] . Moreover, in this phase, as mentioned before recent 
studies indicate that the collective effects could 

be suppressed by matter, which means that the neutrino 
conversion would be given by the standard scenario. This 
would allows us to identify the oscillation effect more eas- 
ily, provided that the number of events for this phase is 
large enough. 

However, it must be stressed that differences or similar- 
ities between 9,, and fluxes must be confirmed (or re- 
futed) by observations. The real SN neutrino data could 
be very different from our expectations! We also should 
keep in mind that the characteristics of SN neutrinos may 
depend strongly on the SN (on the progenitor's mass, 
chemical composition of the core, etc.) (see e.g., [ST]V 

we show the same quantity, R^Nj^^p/N^p), 



In Fig. 



11 



but as a function of the true (input) value of for 
the cases where the true value of {Ey_^) is 18 MeV. A 
similar conclusion can be drawn as before. Unless {Ej^^) is 
significantly different from {E^^ ) , it will not be very easy 
to distinguish p — from p = 0.69. We can also use these 
results to do the opposite. Namely, if we know the mass 
hierarchy by the time of the next galactic SN observation, 
we can try to infer, to some extent, the original values of 
{Ejj^) and {Ei,J) from data. 

we show R{Njy^p/Nt,p) 



In Fig. 
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) as a function of p for 
the cases where the true values of are 18 and 21 

MeV. We see again that unless is significantly dif- 

ferent from {Ep^ ) , it is not easy to identify the oscillation 
effect, especially if the hierarchy is normal and p = 0.69. 
The case of inverted hierarchy, p — 0.0, seems easier to 
be identified. 

we show the ratio R{Np^plNi,p) as a function 
((£:i>J<I>i,J/((E;,,J$^J for the case where 



In Fig 
of LyjLJl 
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the true value of (Ey^ is 18 MeV. This plots indicate 
that if the luminosity of Vf, is larger than that of by 
^ 20 % or so, the cases of normal and inverted mass 
hierarchy can be confused. However, if v^, luminosity is 
significantly larger or smaller (by ^ 50 % or more), then 
it is easier to distinguish the mass hierarchy. On the 
other hand, if the mass hierarchy is known by the time 
the next SN neutrinos are observed, then we can try to 
infer the difference of the luminosities of v„ and v^. 



E. Earth Matter Effect: Shadowing Probabilities 

One potentially interesting possibility is to observe the 
Earth matter effect due to SN neutrinos [551 [75H77] . If 
observed, it could unravel the neutrino mass hierarchy 
and some properties of the SN neutrino fluxes. Suppose 
that SN neutrinos are observed by more than two neu- 
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FIG. 10: The ratio, R{Nii^p/Nvp), as a function of the true 
value of (-E^^) for our reference ANDES neutrino detector. 
The values of the SN parameters not specified in this plot are 
the ones defined in Table Hill 
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FIG. 11: Same as Fig. |10| but as a function of the true value 
of (Ev^). The values of the SN parameters not specified in 
this plot are the ones defined in Table [111] 
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FIG. 12: Same as Fig. [To] but as a function of p for the cases 
where the true value of (-E^^.) = 18 and 21 MeV. 
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FIG. 13: Same as Fig. |10| but as a function of the true value 
of Lp/L^^ = ({Ei^J$pJ/{{E^J^^J for the cases where the 
true value of {E^^ ) = 18 MeV. The values of the SN param- 
eters not specified in this plot are the ones defined in Table 
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trino detectors. If some of them (but not all) receive SN 
neutrinos passing through the Earth interior (shadowed 
by the Earth) but at the same time, there are some other 
group of detectors which receive SN neutrinos without 
passing through the Earth interior (non shadowed by the 
Earth), it would be interesting to compare the results of 
these two group of detectors. 

Currently, several neutrino detectors, which can ob- 
serve neutrinos coming from a galactic SN, are in opera- 
tion. Among them we can highlight Super-Kamiokande 
(SK) 46;, KamLAND g], Borexino IceCube gT] and 
LVD 78 . In this paper we will consider the following lab- 
oratory sites: Kamioka (SK and KamLAND), the South 
Pole (IceCube), Sudbury (SNO-h) [E] and ANDES 7J. 
In Table |V| we show the positions of the detectors we 
consider in this work. 

Following Ref. [52], we have calculated the shadowing 
probabilities for the cases where N {N = 1,2,3,4) de- 
tectors are considered simultaneously. The shadowing 



Site 


Latitude 


Longitude 


Shadowing Prob. 
Mantle (Core) 


Kamioka, Japan 
South Pole 
ANDES 
SNO, Canada 


36.42°N 

90°N 
30.25°S 
46.476°N 


137.3° E 

68.88°W 
81.20°E 


0.559 (0.103) 
0.413 (0.065) 
0.449 (0.067) 
0.571 (0.110) 



TABLE V: Positions of the detectors we consider in this pa- 
per. In the last column, we show the shadowing probability, 
the probability that SN neutrinos will pass only through the 
mantle (indicated as Mantle) or will pass both through the 
mantle and core (indicated as Core). 



probability is defined as the probability that a given de- 
tector or combination of detectors receives neutrinos from 
a galactic SN with Earth matter elfect either by passing 
only through the mantle or both through the core and 
the mantle of the Earth [52, . 

Using the same model of SN distribution in the Milky 
Way considered in Ref. [52], which is based on the neu- 



tron star distribution (reproduced in Eqs. (CI) and (C2) 
of the Appendix C]) we can compute the shadowing prob- 
ability for an arbitrary number of detector positions on 
the Earth. 

The shadowing probability for the most simple case 
where only a single detector is considered is shown in the 
last column of Table |Vl For the case where two detector 
locations, Kamioka and South Pole, are considered simul- 
taneously, we show the results in Table |VI[ We observe 
that the numbers shown in this table agree well with the 
ones found in Table 2 



we 



of Ref. [52]. From Table |\T 
conclude that the probability that only one of these de- 
tectors observes SN neutrinos having passed through the 
Earth is about 72%. 





Earth Matter Effect 




Case 


Kamioka 


South Pole 


Shadowing Prob. 








Mantle (Core) 


(1) 


No 


No 


0.152 (0.832) 


(2) 


Yes 


No 


0.435 (0.104) 


(3) 


No 


Yes 


0.288 (0.065) 


(4) 


Yes 


Yes 


0.125 (0.000) 



TABLE VI: Earth shadowing probability for the case where 
two detectors at Kamioka and South Pole are considered. 



Next in Table [VTl] we show the Earth shadowing prob- 
abilities for the case where we consider three detectors: 
at Kamioka, South Pole and ANDES sites. From this ta- 
ble, we see that the probability of having at least one of 
these detectors observing SN neutrinos passing through 
the Earth while at least one of the other two sees them 
non shadowed by the Earth is 96%, which is about 30% 
larger than the case above with two detectors, one at 
Kamioka and the other at the South Pole. One can also 
compare our three detector combination with any two de- 
tector combination found in Table 2 of Ref. 52 where the 
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Earth Matter Effect 




Case 


Kamioka 


South Pole 


ANDES 


Shadowing Prob. 
Mantle (Core) 


(1) 


No 


No 


No 


0.024 (0.767) 


(2) 


Yes 


No 


No 


0.388 (0.105) 


(3) 


No 


Yes 


No 


0.034 (0.061) 


(4) 


No 


No 


Yes 


0.128 (0.063) 


(5) 


Yes 


Yes 


No 


0.106 (0.000) 


(6) 


No 


Yes 


Yes 


0.254 (0.003) 


(7) 


Yes 


No 


Yes 


0.047 (0.000) 


(8) 


Yes 


Yes 


Yes 


0.020 (0.000) 



TABLE VII: Earth shadowing probability for the case where 
detectors at Kamioka, South Pole and Andes sites are consid- 
ered. 



by comparing the yields of two (or more) detectors if 
only some (not all) of them receive SN neutrinos passing 
through the Earth's interior. See [S31 [75H77] for the 
previous works. 

In this paper, we focus on the possible Earth matter 
effect for Vf. due to the larger number of expected events. 
Since we know now that ^13 is not too small, sin^ 2^13 ~ 
0.09-0.1 [nHHIMlISi, it is expected that Earth matter 
effect can only be large for the normal mass hierarchy in 
the standard scenario we consider in this work [53j . So 
we will assume here only the normal hierarchy in such 
a standard scenario. Note that for the case of Vf. in the 
normal mass hierarchy, collective effects, shock wave, etc. 
are in general expected to be small. 



largest probability for having one detector shadowed and 
one non shadowed is 87.2%, occurring for the Pyhasalmi 
and South Pole sites. 

In Table IVIIII we show the case where four detectors 
at Kamioka, South Pole, ANDES and Sudbury sites are 
considered. With four detectors, the probability that at 
least one of the detectors have the Earth effect and at 
least one of the others have not, increases to 98%. We 
also found that the probability that at least one of these 
sites receives a SN neutrino flux which passes the core of 
the Earth is not very small, ^ 34%. 

F. Quantifying the Earth matter effect: Comparing 
the detectors with and without Earth matter effect 



If SN neutrinos reach the detector after passing 
through the Earth matter, assuming the impact of non- 
zero 013 is not very significant for the Earth matter effect 
itself (though the impact of ^13 is expected to be large 
for the oscillations inside the SN, affecting significantly 
especially for the case of the inverted mass hierar- 
chy), the SN flux spectrum is expected to get modified 
as follows [53] . 



Let us now try to quantify the Earth matter effect 
which could be observed in a model independent way where 



f{E) 



\U, 



e2| 



It/. 



(8) 



where Vek (fc = 1, 2) are the elements of neutrino mixing 
matrix which relate flavor and mass eigenstates and 

=P®(;>fe^^^„L), (A; = 1,2), (9) 

are the probabilities that a mass eigenstate (fc — 1, 2) 



entering the Earth will be detected as v^, at the detector, 
after traveling the distance L inside the Earth. 

If we take the difference of SN spectra given in Eqs. 
([7]) and Q with and without Earth matter effect. 



1 



cos 2Q^ 



mE) - 1](.L -P?e) {fUE) - FliE)} 



ip?e-cl,){FliE)-FliE)} 



(10) 
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Case 



(1) 



(2) 
(3) 
(4) 
(5) 



(6) 
(7) 
(8) 
(9) 
(10) 

(11) 



(12) 
(13) 
(14) 
(15) 



(16) 



Earth Matter Effect 



Kamioka South Pole ANDES 



No 



Yes 
No 
No 
No 



Yes 
Yes 
Yes 
No 
No 
No 



No 
Yes 
Yes 
Yes 



Yes 



No 



No 
Yes 
No 
No 



Yes 
No 
No 
Yes 
Yes 
No 



Yes 
No 
Yes 
Yes 



Yes 



No 



No 
No 
Yes 
No 



No 
Yes 
No 
Yes 
No 
Yes 



Yes 
Yes 
No 
Yes 



Yes 



Sudbury 



No 



No 
No 
No 
Yes 



No 
No 
Yes 
No 
Yes 
Yes 



Yes 
Yes 
Yes 
No 



Yes 



Shadowing Prob. 
Mantle (Core) 



0.008 



0.206 
0.034 
0.001 
0.016 



0.205 
0.000 
0.282 
0.163 
0.000 
0.127 



0.091 
0.047 
0.011 
0.012 



0.008 



0.657) 



0.105) 
0.061) 
0.063) 
0.111) 



0.000) 
0.000) 
0.000) 
0.003) 
0.000) 
0.000) 



0.000) 
0.000) 
0.000) 
0.000) 



0.000) 



TABLE VIII: Earth shadowing probability for the case where detectors at Kamioka, South Pole, ANDES and Sudbury sites 
are considered. 
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FIG. 14: Neutrino Oscillogram, iso-contour plot of the ratio P®(i^i — >■ «^e)/ci2 in the plane of the neutrino energy and the nadir 
angle. Note that unity for this ratio corresponds to the case where the Earth matter effect is absent. 



where p{E) sa c^2 "^^^ assumed to get the last expression. 
We also observe that the term proportional to |C/e3p can 
be dropped because it will only contribute to about 7% 



of the first term. As one can see from Eq. (10 1, in order 



to observe the Earth matter effect, the deviation of 
from must be large enough and at the same time, 
the difference between the original spectra of (E) and 
(E) must be also significantly large. 
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FIG. 15: The fractional difference of the expected neu- 
trino flux spectra with and without the Earth matter ef- 
fect, AFiiJFi,^ where AFi,^ = F^^{E) - Fi,^{E), given in Eq. 
( |10[ |. The SN neutrino parameters are the same as assumed 
in Fig. [51 



In order to have some idea about the magnitude of the 
Earth matter effect as functions of the neutrino energy 
and the incident angle of SN neutrinos, we show in Fig.[T4| 
the iso-contours of the quantity P®{i>i — >■ i'e)/cl2 in the 
plane of the nadir angle, 6'nadir, and the neutrino energy 
where P®(Pi — > i/g) was obtained by numerically inte- 
grating the neutrino evolution equation using the matter 
density profile of the Earth based on PREM (Preliminary 
Reference Earth Model) [79] . Note that any deviation of 
P®(i/i — >■ i'f,)/ci2 from unity implies the presence of the 
Earth matter effect. 

This plot is similar to the so called "neutrino oscil- 
logram" studied in detail in Ref. [80] in the context of 
atmospheric neutrinos. As we can see from this plot, the 
Earth effect is as expected strongest when neutrinos pass 
through the Earth's core (corresponding to ^nadir ^ 33°) 
and for higher neutrino energies. Here 6'nadir is defined 
such that 6'nadir = 0° corrcspouds to the case where SN 
neutrinos arrive at the detector from the other side of the 
Earth passing the center of the Earth and 6'nadir = 90° 
corresponds to the case where neutrinos come from the 
horizontal direction. 

In Fig. [15] we show the fractional difference of the 
flux spectra with and without the Earth effect defined 
by AFp^/Fp^ where AFp^ is given by Eq. (10) for for 



three different path-lengths in the Earth, L = 1000 km 
(solid blue curve), L = 5000 km (dotted red curve) and 
L = 10000 km (dashed green curve). From this figure, 
we can see that strong Earth matter effect is expected in 
the higher energy range. We must note, however, that 



as the energy becomes higher, the number of events gets 
smaller, so that in order to identify the Earth matter ef- 
fect, both the Earth effect (difference of probabilities with 
and without matter effect) and the number of events in 
the relevant energy range must be large enough. 

Since we cannot compare the number of events at the 
same detector with and without Earth matter effect, we 
need two or more detectors to be able to conclude some- 
thing on the matter effect. For simplicity and for the 
sake of discussion, let us consider only two detectors at 
two different sites, say, one at Kamioka (SK) and other 
at ANDES. 

Suppose that the arrival of the SN neutrinos at the 
ANDES detector is shadowed by the Earth, while at SK 
detector at Kamioka is not. Then, to some extent, within 
the statistical and systematic errors, one can try to infer 
the expected SN spectra at SK from the observed ones 
at ANDES (or vice versa) as if the SK detector were 
also shadowed like ANDES. If both detectors receives SN 
neutrinos without Earth matter effect these two spectra 
must agree with each other but with the matter effect, 
they are not expected to coincide exactly. 

In order to see the presence of Earth matter effect, the 
combination of SK and ANDES must be able to distin- 
guish, for a given SN model, vacuum from matter event 
distribution. To illustrate that, we show in Table [IX] the 
expected number of events in SK for four different en- 
ergy bins for vacuum and matter with L = 1000 km, 
for a variety of SN parameters. If we assume that the 
shadowed SN neutrino spectra distribution of events at 
SK can be provided by ANDES with an uncertainty in 
each bin equal to the statistical uncertainty in ANDES 
we can estimate in which of these cases one can establish 
the presence of matter effects. In fact, since the first bin 
-E < 30 MeV is not sensitive to matter effects, one can 
use it as a control bin to evaluate the accuracy of the 
vacuum distribution provided by ANDES data. For this 
study, we consider SN event at 5 kpc from the Earth as 
we need larger number of events for ANDES. 

To give a quantitative idea of the discrimination power, 
we have assumed, for each SN model in Table |IX| the 
vacuum distribution to be the hypothesis to be tested 
with data at 1000 km. We use the last two energy bins in 
order to estimate the deviation from vacuum in terms of 
standard deviations. These numbers are given in the last 
column of Table [TX] For most cases, matter effects can be 
established in more than 2 a. We have also verified that 
if < 18 MeV we cannot distinguish matter effects 

for any value of the other SN parameters considered in 
Table [!xl at least if L^-JL^^ = 1. 



IV. DISCUSSIONS AND CONCLUSIONS 

The ANDES laboratory, if constructed, will be the 
first deep underground laboratory in the Southern Hemi- 
sphere. It can offer the possibility to explore rare physics 
events that can profit from the low natural background 
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£; < 30 MeV 


30 < E/MeV < 40 


40 < E/MeV < 50 


£ > 50 MeV 


Case 


Compatibility 


Vacuum (observed) 
1000 km (prediction) 


18159 ± 135 
18132 ± 374 


4973 ± 71 
5065 ± 198 


2032 ± 45 
1908 ± 121 


889 ± 30 
700 ± 74 


(B.J = 22 MeV 

= /3e = 4 


3.1a 


Vacuum (observed) 
1000 km (prediction) 


17395 ± 132 
17370 ± 367 


5785 ± 76 
5858 ± 213 


2583 ± 51 
2483 ± 139 


1147 ± 34 
988 ± 87 


(B.J = 22 MeV 

/3, = 4 /?e = 3 


2.2a 


Vacuum (observed) 
1000 km (prediction) 


16031 ± 127 
16011 ± 352 


6674 ± 82 
6728 ± 228 


3594 ± 60 
3541 ± 166 


1978 ± 45 
1917 ± 122 


(E,,,) = 22 McV 
/3, = 4 /?e = 2 


0.7a 


Vacuum (observed) 
1000 km (prediction) 


16863 ± 130 
16837 ± 361 


5722 ± 76 
5787 ± 212 


2864 ± 54 
2731 ± 145 


1604 ± 40 
1321 ± 101 


{E^J = 22 MeV 


3.2a 


Vacuum (observed) 
1000 km (prediction) 


15499 ± 125 
15479 ± 346 


6611 ± 81 
6657 ± 227 


3875 ± 62 
3789 ± 171 


2434 ± 49 
2250 ± 132 


{E^J = 22 MeV 

= 3 /?e = 2 


1.7a 


Vacuum (observed) 
1000 km (prediction) 


14790±122 
14766±338 


6388± 80 
6419±223 


4089± 64 
3971± 175 


3059± 55 
2701± 145 


= 22 MeV 


2.8a 


Vacuum (observed) 
1000 km (prediction) 


17686±133 
17655±370 


5240±72 
5343±203 


2439±49 
2272±133 


1285 ±36 
990±88 


{E,J = 24 MeV 
= /3e = 4 


4.3a 


Vacuum (observed) 
1000 km (prediction) 


16922±130 
16892±362 


6052±78 
6136±218 


2990±55 
2847± 148 


1543±39 
1278±100 


(E^^) = 24 MeV 
= 4 /3e = 3 


3.1a 


Vacuum (observed) 
1000 km (prediction) 


15557±125 
15533±347 


6941±83 
7006±233 


4001±63 
3905±174 


2374±49 
2207± 131 


(E,^) = 24 MeV 
/3, = 4 /?e = 2 


1.7a 


Vacuum(observed) 
1000 km (prediction) 


16441±128 
16409 ±356 


5858±77 
5928±214 


3174±56 
3007± 153 


2022±45 
1625±112 


(B.J = 24 MeV 
= = 3 


4.0a 


Vacuum (observed) 
1000 km (prediction) 


15077±123 
15051±341 


6746±82 
6797±229 


4185±65 
4065±177 


2853±53 
2554±141 


(B^J = 24 MeV 
= 3 /3, = 2 


2.5a 


Vacuum (observed) 
1000 km (prediction) 


14439±120 
14410±334 


6400±80 
6430±223 


4248±65 
4116±179 


3402±58 
2948±151 


(B^J = 24 MeV 
= /3. = 2 


3.3a 



TABLE IX: Number of inverse /3-decay events expected at the SK detector (1.7 xlO^^ free protons) for a SN happening at 5 
kpc from the Earth for £; < 30 MeV, 30 < E/MeV < 40, 40 < E/MeV < 50 and E > 50 MeV for the case of vacuum and 
matter effect with a basehne of 1000 km and various SN parameters. We assume i.- /L.^ — 1 and (Ej^^) = 15 MeV. We assume 
that the SK detector receives SN neutrinos without the Earth matter effect whereas the ANDES neutrino detector receives 
them after travehng 1000 km inside the Earth. The numbers in the row indicated as vacuum are the ones to be observed at 
SK detector which must be compared with the theoretical prediction at SK inferred from the observed number of events at 
ANDES neutrino detector. In the leist column we point out in how many a the 1000 km observation is distinguishable from 
vacuum. 



environment because of the overburden of ~ 4.5 kmwe. 
In particular, dedicated neutrino physics and dark mat- 
ter search programs, that could also benefit from the its 
unique geographic location, could be envisaged. See ^ 
for updates on the status of the laboratory. 

In this work, we have studied the potential of a few kt 
liquid scintillator neutrino detector at the ANDES un- 
derground laboratory for neutrino astrophysics and geo- 
physics. Since there are very few nuclear reactors in 
South America, the location of the ANDES laboratory 
is specially suitable for geoneutrino observations. More- 
over, due to thick continental crust around the labora- 
tory, higher geoneutrino event rate is expected, substan- 
tially larger than at Kamioka and Gran Sasso, which by 
itself would be interesting to confirm experimentally. 

Concerning the observations of neutrinos coming from 
SN, the ANDES neutrino detector could play an impor- 
tant role. First of all, because of the small event rate of 
the nearby SN (within a few 10 kpc from the Earth or 



so) , having as many large neutrino detectors as possible is 
highly desirable. Furthermore, because of the location, 
the presence of the ANDES detector will increase the 
chance of observing the Earth matter effect by combin- 
ing the signal at the ANDES with other detectors in the 
Northern Hemisphere. The ANDES neutrino detector 
could also integrate the international SN watch network 
SNEWS [27]. 

We have focused here on galactic SN neutrinos and 
geoneutrinos. However, with such a liquid scintillator 
detector one could also try to study solar neutrinos in a 
similar fashion as Borexino [5] and KamLAND [ST] have 
done and as SNO+ [15] plans to do. In fact, Borexino 
has shown how pep, CNO and perhaps even pp solar 
neutrinos can be accessible by such a detector. There is 
also an interesting proposal to search for indirect dark 
matter signals through neutrinos which are coming from 
the dark matter annihilation in the Sun in these type of 
detectors (see e.g. ^82]). 
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energy, Cy = 0.04 and Ca = 1.27/2. To take into account 
the loss of proton energy, we calculate the quenched pro- 
ton energy, T', by 



dT 



keidT/dx)' 



(B2) 



where kg is called Birks constant. The numerical values 
of (dT/dx) which describe the energy loss of proton in 
scintillator were taken from Ref . [86] . 

The event number distribution dN/dT' is calculated 

by 



dN 
dT' 



°° dF da 

dE , B3) 



where i?min is the minimum energy of neutrino which can 
produce recoil proton with the kinetic energy T, given by 



Appendix A: Calculations of the number of event 
induced by the inverse beta decay reaction 

The number of event induced by the inverse beta decay 
reaction j/g + p — > ri + e+ , is given by 



nOO 

p dE F^^{E)a^^p{E) (Al) 



where £^min = 1-806 MeV is the threshold of this reac- 
tion, Np is the number of free protons in the detector and 
E is the observed energy. For simplicity, we assume per- 
fect energy resolution, which is a good approximation for 
the results presented in this work. While the accurate 
absorption cross section of S'e on proton, ajj^p{E), can 
be found in Ref. [83!, we only quote the simple approxi- 
mate analytic expression given in this reference, which is 
sufficient for our purpose, 



PeEeE 

xlO"''^ 



-0.07056+0.02018 In £-0.001953 In^ E 
21 



cm 



(A2) 



where E,^ — E^^nin^mp) ~ £'—1.293 MeV (r7i„, are 
mass of neutron and proton, respectively) , and Pe is the 
momentum of positron, all the energies should be given 
in McV. 



Appendix B: Calculations for proton- neutrino 
scattering 

The differential cross section for v + p ^ t^+P, for any 
neutrino flavor, is given by [68l [84] . 



da 
dT 



GpTUp 



1 



mpT , 2 



1 



2£2 



(Bl) 

where Gp is the Fermi constant, rup is the proton mass, 
T is the kinetic energy of recoil proton, E is the neutrino 



1 



T + JT{T + 2mp 



nipT 
2 



(B4) 



We further taken into account the detector energy resolu- 
tion to convert dN/ dT' in Eq. ( B3 ) into the ones actually 
observed, dN/dTquc as 



dN 
dT 



que 



pOO 1 TV T- 

dT'—R{T^^,X), (B5) 



where R{Tquc-,T') is the resolution function given by 



RiTquc,T') — 



1 



27r(7 



cxp 



2a2 



(B6) 



with the resolution assumed to be 5 %/ y^T' /MeV [69] , 



Appendix C: Some details of the calculations of 
shadowing probabilities 

For definiteness, we use the same galactic supernova 
distribution model considered in Eqs. (1) and (2) of 
Ref. [S2], which is based on the expected distribution 
of neutron stars in the Galaxy ^85] given by 



(Jsn{i") oc exp 



1.25 kpc 



(CI) 



where (Tsn('") is the surface density of the core-collapse 
SN events as a function of the radial distance from the 
galactic center. We also take the same vertical distribu- 
tion of the SN events as 



RsNiz) oc 0.79 exp 
+0.21 exp 



\212 pc 



636 pc 



(C2) 



18 



0.10 
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Pulsars 
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NS: <d> = 10.7 kpc 
Pulsar: <d> = 10.3 kpc 
laSNe: <d> = 12.2 kpc 




0.00 
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D (distance to SN from the Earth) [kpc] 

FIG. 16: Probability of SN distribution in the Milky Way 
based on the distributions considered in Ref. [55]. The red 
solid, dashed blue and dotted green curves indicate, respec- 
tively, neutron star, pulsar and type la SN. In this work we 
consider the neutron star distribution (red curve) as our ref- 
erence SN distribution. 



where z is the vertical distance from the galactic plane 
so that the SN distribution ns-t^{r,z) oc crsN(?')-RsN(-z)- 



In Fig. [16] we show the expected SN probability distri- 
bution in the Milky Way based on these distributions by 
the solid red curve. For comparison, we show the other 
distribution considered in Ref. [52j which are based on the 
pulsar (blue dashed curve) and type la SN (green dotted 
curve). We note that in calculating the shadowing prob- 
abilities rotation of the Earth was taken into account in 
the same way as in [521 . 
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